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INTRODUCTION 


Eddy current inspection of graphite epoxy filament wound materials has been 
investigated by several researchers 1 ' 12 . The results of these investigations show 
considerable promise for eddy current testing to locate flaws such as cut tows and 
impact damage in critical structures fabricated from graphite filaments. Detection of 
delaminations in such materials is more difficult in that voids in the conducting media, 
such as observed in kissing debonds, do not always exist naturally. Loaded samples may 
be more easily inspected however, because of possible displacement of the laminae 
surfaces to create a non-conducting volume that is observable using eddy currents. 

Earlier work of this research group resulted in the developement of a robotic scanning 
eddy current inspection cell which is located in Building 4702 at Marshall Space Right 
Center. A Macintosh FX computer with a Symbolics Mclvory board set provided the 
computational and control function. However, eddy current sensing was performed 
with either a Smarteddy (PC platform) or a HP 4193A Impedance Meter. Data 
conversions to get the data into the Macintosh were either cumbersome in the former 
case or extremely slow in the latter case. Consequently this work concentrated on 
designing and building an eddy current interface for the Maclvory which could replace 
both of the above instruments for acquiring eddy current signals. 

The basic concept for this research activity is presented in Figure 1. Manual inspection, 
particularly for either large filament wound items or large volumes of small cylindrical 
items, would be a very inefficient and tedious process. Eddy current waveforms of the 
composite filament componentsare rather noisy since the filaments themselves provide 
a time varying signal as the probe scans over them. Any other motions caused by 
manual scanning manipulations would then add to the already modulated signal 
normally expected and would make interpretation more difficult. Consequently this 
research has concentrated on automated systems approach. Adding an expert system to 
the interpretation is still a major goal once a correct set of interpretation procedures have 
been accomplished. 



Figure 1. Overall Systems Concept for the Robotic Eddy Current Inspection Workcell. 
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EXPERIMENTAL 
Workcell Components 

The Robotic Scanning Laboratory at the NDE Laboratory at MSFC consists of the 
workcell items as shown in Figure 1. An Indellidex 550 robot is used to orient various 
eddy current probes onto the part under test and perform scaning operations. The five 
degrees of freedom available in the robot arm are augmented by a DCI turntable which 
rotates the item under test. All control functions are programmed through the Macintosh 
system running Labview. A sketch showing how the eddy current transducer atached to 
the robot end-effector approaches the item under test is shown in Figure 2. Due to the 
flexibility of the robotic end-effector a variety of probes and trajectories are available for 
NDE applications. 


2 




Figure 2. Position of Robot Arm in Approaching Part Under Test 



Notice that the transducer shown in Figure 2 is an E-probe. As part of this work, we 
have fabricated several types of probe assemblies to be used with the robot system, 
including a horseshoe, an E-probe, a cup probe and a pancake probe, the dimensions of 
these probes are shown in Figure 3. Each probe uses ferrite material to concentrate the 
magnetic flux into the part under test. The ferrite material for all probes except the 
pancake probe was a high permeability ferrite component from Magnetics, Inc. The 
ferrite rod used to prepare the pancake probe was obtained through a local electronics 
parts supplier. 
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Figure 3. Various Eddy Current Probes Used to Inspect Graphite Epoxy Components 
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Eddy Current Instrumentation Concepts 

The earlier work used the HP 4 193 A to generate and acquire the eddy current 
signals for the graphite epoxy samples. Since the maximum data sampling rate was very 
slow (7 readings per second), a major effort was made to develop the capability for 
faster data sampling in a mode that was compatible with the Macintosh computer 
system. In order to accomplish this task, we chose to fabricate our own circuitry for 
acquiring the eddy current signals and thus to use the HP 4193 only to generate the 
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alternating current needed to generate eddy currents. Figure 4 shows the instrumental 
concept employed for both both generating, decoding phase and magnitude, and 
acquiring the eddy current signals. 


Figure 4. Schematic of Eddy Current Source and Detection System Used in this Work. 



The circuitry used here is very similar to other systems and uses the XAR 2208 RMS- 
DC converter devices to generate a dc level signal from the ac signals observed during 
the eddy current scan. Using two XR 2208's allows both the phase and the magnitude to 
be digitized and acquired with a GW Instruments A/D converter. The GW Instruments 
A/D uses the SCSI bus to transfer the digitized data to the Macintosh. 
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Figure 5. Schematic of Bridge Circuit used to acquire Eddy Current Signals. 



Test Specimens 

Cylindrical samples of graphite epoxy with two different thicknesses were used 
for these studies. These samples were prepared by MSFC personnel in their filament 
winding laboratory in Building 4707. The first specimen was a four ply, 12 inch in 
diameter cylinder, representative of thin specimens. The second specimen was 24 plies 
thick with the same diameter. The layout of the manufactured defects for the two 
specimens are shown in Figures 6 and 7 on the next page. The depth of penetration for 
these specimens is shown in Figure 8. This plot allows one to determine sensitivity for 
flaw detection according to ply number in the composite. 
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Figure 6. Manufactured defects for four ply test cylinder. 
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Figure 7. Manufactured Defects for Thick specimen 








Figure 8. Depth of Penetration of Graphite Epoxy at 3.2 Mhz 
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Quantitative Modeling Of Responses In Eddy Current Probes 

The ability to characterize unknown defects detected by the eddy current method 
is limited by the lack of adequate theoretical models. Basically, such models should take 
the place of very difficult and expensive experiments in developing the knowledge base 
necessary to approach the inverse problem that given a particular eddy current response, 
characterize the defect. Extensive developments of the finite element method for eddy 
current problems has been very successful for two-dimensional or axisymmetric 
geometries. However, the difficulties occur in the FEM three-dimensional case due to 
the large number of nodes required and the necessity for three non-zero components of 
the magnetic vector potential. 

The FEM model should be most efficient for rather large-size defects because the 
number of nodes necessary to map the fields adequately is then relatively small. The 
model developed by Burrows 33 gives the response to defects of spheroidal shape with 
dimensions small compared with the skin depth in the materials. The boundary-element 
integral equation model of Hower 34 has allowed consideration of defects of various 
shapes but with the same size restriction. Still needed is an efficient method for 
calculating the effect of three-dimensional defects having dimensions up to one to two 
skin depths. 

The model developed in this work can be considered as a hybrid approach 
utilizing several available techniques. The method involves separation of the problem 
into three parts: (1) Use a commercially available finite element software, Ansoft 
Maxwell in our case, to calculate the incident fields of the horseshoe probe, but in the 
absence of any defect. (2) Compute the total fields on the surface of the defect. (3) 
Determine the resulting change in probe coil impedance from both unperturbed fields 
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and equivalent source on the defect surface using the Lorentz reciprocity theorem. The 
mathematical foundation for this operation follows below. 

The change in coil impedance AZ due to the flaw can be found in terms of an 
integral of the field quantities E and Eq (the fields at the position of flaw with and 
without the flaw being present). These formulas have been derived by Auld 29 and 
Zaman, et al 30 . 

Let us assume that within a linear, isotropic medium, which is not necessarily 
homogeneous, there exist two sets of sources Jj, Mj that are allowed to radiate 
simultaneously or individually inside the same medium at the same frequency and 
produce fields Ej, Hj and E 2 , H 2 respectively. The field equations are then: 


V x H j = J j + jcoeEj (1) 

V x Ej = -Mj + jcoeHj (2) 

VxH 2 = J 2 +jcoeE 2 (3) 

V x E 2 = -Mj + j(oeH 2 (4) 


We multiply Equation (1) scalarly by E 2 and equation (4) by Hj and add the resulting 
equations. Similarly Equations (2) and (3) to give: 

V ( E 2 x Hj) = + jco£EjE 2 + jcoeHx^ + E 2 JJ + HjM 2 

V ( Ej x H 2 ) = + jco£EjE 2 + jco8HjH 2 + Ej J 2 + H 2 Mj 

where the left-hand term has been simplified by the vector identities. A subtraction of the 
former equation from the latter yields: 

V ( Ej x H 2 - E 2 x Hj ) = E 2 Jj - E t J 2 +H!M 2 - 

At any point for which the fields are source-free (J=M=0) this reduces to: 

V ( Ej x H 2 - E 2 x Hj ) = 0 

which is called the special case of Lorentz reciprocity theorem. Its integral form is: 


J (E 1 xH 2 -E 2 xH 1 )dS (5) 

s 

Expression for the change in coil impedance (AZ) are obtained by application of above 
reciprocity theorem. 

The simplest geometry appropriate to this problem is the case shown in Figure 9. 
It consists of a conductive body B, a defect and a eddy current horseshoe probe. Identify 
the multiply connected region bounded by the surface Sj of the defect, the closed surface 

S s surrounding the source, the cable cross-section S c and a boundary at infinity S„„. Let 
the volume enclosed by this surface be V and Eq, Hq denote the time harmonic fields 
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when the defect dose not exist and E, H denote the field actually existing with the defect 
present 



Since the volume V contains no source currents (the scattered fields are 
considered as being due to surface currents on S^) and the material within V is identical 
with or without the presence of defect, the Lorentz reciprocity theorem is valid, 


J (E 0 xH-ExH 0 )dS = 0 


where dS is an element of surface which is normal to and directed into the bounded 

volume. The surface S is the sum of surfaces which include S„,, S c , S s and Sj. 
However, 


J-M+M-o 

S S, s, s d 

Since the integrand goes to zero on S„, as this partial surface is allowed to move to 
infinity, and it also goes to zero on the source enclosure S s , if it is assumed to be 
perfectly conducting. There remains the feed line cross-section S c and the surface S<j 
enclosing the defect. Hence, we have 


J (EoxH-ExH 0 )dS= J (EoxH-ExH 0 )dS 


( 6 ) 


10 




Assuming a coaxial cable feed line, the surface S c appears as shown in Figure 10. If the 
field quantities in Equation 5 are defined as those produced (either in the absence or the 
presence of the defect) by a drive current I at the surface S c , then the magnetic fields on 
S c are 


in terms of the feed line current /; and the feed line voltages in the absence and presence 
of the flaw are given in terms of the electric fields by 

J Ej-dr = -ZI, | Eq,. dr =- Zq 

where Z designates the circuit impedances at this point in the line. Note that the field 
distributions, which are those characteristic of a coaxial line, are unchanged by the 
presence of the defect because of the way these fields are normalized, only the electric 
field is changed in magnitude. From this, it is found 


J ( E 0 x H - E x H 0 ) dS = J (E 0 xH <p -E r xH 0<p )dS = l2(Z0-Z) (7) 



— =► H-field 
► E-field 

Figure 10. Distribution of electrical and magnetic fields inside a coaxial transmission 
line. 

Substitute Equation (7) into Equation (6), gives 
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( 8 ) 


i 2 (Zo - z) = - J 


s. 


(EoxH-ExH 0 )dS 



(9) 


Using Gauss's theorem, 

AZ = -I J [H (V x E 0 )- E 0 (V x H) - H 0 (V x E) + E (V-x H 0 )] dV 
1 K i 

AZ = -jrJ [H (-j®Mo)- Eq (jcoeE) - H 0 (-jtotloH) + E -(jtteE 0 )]dV 

1 Yt 

AZ = j Eo . Edv 

1 K 


But 


£ = £ — — 
c c 0 to ’ 


* £ — £ = — — 


XZ = ^jE Q *EdV (10) 

1 v j 

where E and Eo are the electric fields with and without the defect present. Equation (10) 
is the formula that we use in this work. 

The case of a defect whose permittivity is near that of the host material may be 
examined in greater detail by allowing the field in the presence of the defect E to be 
approximated by the unperturbed field Eq over the volume of the defect. For defects 
small enough that the fields do not vary greatly over its volume, the expression may be 
further approximated by using just the value of the field at the position of the centroid of 
the defect. 


AZ = 



( 11 ) 


Similar derivations for the change in impedance due to a flaw have also been 
developed by Bahr 31 and by Auld ^ 2 . The expression derived here in Equation (11) 
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corresponds to a low frequency approximation, and the fields are already known for the 
unflawed case. It can provide a simpler first order approximation for the change in 
complex impedance. The approximation that E = Eq is more accurate for case in which 
the e and s of the defect differ only slightly from that of host material. In this research 
work, the FEM model can calculate both the E and Eq, so Equation (10) is used to solve 
for AZ. 


Traditionally, eddy current nondestructive evaluation methods have been used 
only to provide qualitative information regarding the presence or absence of a defect. 
However, the measured data also contain quantitative information about the defect size. 
From the finite element method and electromagnetic reciprocity relation, the impedance 
change at the probe terminals is formulated in terms of an electromagnetic field integral 
over a surface surrounding the defect. In this section the two-dimensional finite element 
model of E and horseshoe probe on the single layer isotropic material and also the 
calculation of the fields inside the materials is presented. 

Figure 1 1 shows the horseshoe probe as it passes over a defect. Although the 
two-dimensional case is not a realistic model of any actual defect configuration, it does 
serve to bring out the general nature of the probe-defect interaction problems. 



X-axis 

Figure 1 1 Eddy current horseshoe probe scan through the specimen with a defect. 
The defect response is calculated from the value of the magnetic vector potential. 
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Z = f$Adt 


where I s is the source current. Because the magnetic vector potential has only one 
component in the direction of the current for one turn of a coil have a radius r, the 
impedance is: 



ja)2jtr i A i 

T s 


where A[ is the magnetic vector potential at rp Integrating Equation (12) over an 
elemental area with N s turns per unit area the expression becomes: 

7 j co 2 tw c - A C i N s A 

~ Is 


where r c j is the distance from axis to the centroid of the element, A c j is the magnetic 
vector potential at the centroid of the element and Dj is the area of the element 
respectively. Summing Equation (13) over the elements in the coil's cross section and 
noting that Js=N s I s , the coil impedance is: 


z 


_ jcoln] s 
coil i 2 

1 S 


X r cAAi 


i=l 


(14) 


In three-dimensional geometries, this method can not be used because it assumes 
that the magnetic vector potential is constant along the circumference of the coil. 
However, the Lorentz reciprocity theorem allows us to calculate the three-dimensional 
defect response of the horseshoe probe (uniform field probe). 

Electromagnetic fields in a uniaxial conductor have been studied in order to 
evaluate the eddy current testing of composite material such as graphite epoxy. It is 
important to understand how eddy currents behave in anisotropic material, especially 
uniaxial conductors (graphite fibers) in which the conductivity in the axial direction may 
be several hundred times greater than the transverse conductivity. Ansoft Maxwell 
Solver can model the anisotropic material if we set up the material parameters as a 
tensor or a function of position and frequency. But in this research work, all the 
parameters of graphite epoxy are unknowns. 

Theoretical models of eddy current cylindrical air core probe on graphite epoxy 
have been developed by Bowler 22 and Beissner 33 . These solutions are given in terms 
of half-space scalar potential and dyadic Green’s functions. Following Bowler, the 
corresponding electric and magnetic fields are: 
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where 


E. = Jw. [' V x yX + V x V x JX + (1 - - -j)y, ] 

o. 


H. 


-v^-y.vx+l {vxyX 
dy n $, 


V 


t 


~ d „d 

X„ V 2 

dx n dz 


Y n is transverse electric (TE) potential, Y n " is the transverse magnetic (TM) potential 
for layer n. The potentials satisfy the equations 


and 


v 2i K 4 


li T 
( 5„‘) 2 ” 


0 





= 0 


where TE and TM skin depths are: 





The standard approach to the calculation of the impedance of an eddy current 
probe is through Auld’s 29 reciprocity integral. In the composite material, the integral 
involves certain vector products of the probe field in free space and the fields on the top 
(z=0) and bottom surface (z=-t) of the layered structure, when the TE and TM scalar 
potentials are introduced, the TM component does not contribute to the impedance 


15 



integral. Thus we need only the transmitted and reflected TE potential to calculate the 
probe impedance change. According to Beissner [40], the integral leads to 


AZ = 


2jco r z**x 

wl— 


a(*,o) 


)z • k xa(-k, o) - r(fc)z • kxa s (-it,- 


where a s (k,z) is the Fourier transform of the vector potential, R and T are the reflection 
and transmission coefficients, which defined by 

where 

u,(k,z) = — —rrz»k xd\k,z) 

Hok 

and up(k,0) is the Fourier transform of the reflected TE potential at Z=0, u-p(k,-t) is the 
transmitted field at the bottom of the specimen (z=-t). 

This model is equally applicable to any material that can be treated as a 
multilayer anisotropic conductor. The calculation procedure of Equation (14) can be 
modified for the different shape probes in the feature work. 

ANSOFT MODELS 

During this phase of the eddy current reserach, more procedures for utilizing 
Ansoft Maxwell finite element software were developed. The following simulation was 
developed for a 4-ply graphite epoxy filament wound component using the following 
values: 


Graphite conductivity - 14, 300 ohms^m* 1 

Ferrite conductivity - 1 ohms^nr 1 

Ferrite permeability - 2, 500 

Ply thickness 0.052 inches 

Ferrite dimensions - see chart below 

Flaw dimension - 0.1 x 0.052 inches 2 

Frequency - 3.2 megahertz (optimized for cuto tow defect) 
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Figure 9. Geometry of finite element model for internal flaw (cut tow) 
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Figure 10. Finite element model showing magnetic field lines during eddy current scan 

when probe is away from the flaw. 







Figure 11. Magnetic flux lines when probe is over part. 



The Ansoft software provides the capability to generate a macro which allows 
one to simulate a scanning motion across the part being inspected. During each step 
across the part one can embed into the macro particular calcuations such as total 
impedance at the coil, J or M values along a line or within a region. For example in 
simulating the above scan, the information normally observed is the impedance response 
of the probe. By including a macro to increment in the x direction 0.05 inches each step 
within the mesh generating program and then including within the eddy current solver 
program the calculation for the X R and X L one obtains the following charts: 

Figure 12(a). Chart showing expected eddy current resistance measurement when 
scanning part shown in Figures 10-11. 
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Figure 12(b). Chart showing expected eddy current reactance measurement when 
scanning part shown in Figures 11-12 


Inches Scanned 


This simulated response can be compared with actual data obtained from an eddy 
current inspection.The experimental data obtained from the robotic eddy current 
workcell provides a much noiser signal. The following waveforms are obtained from 
scanning the thin filament wound case over regions of no known defects and regions 
with known defects. 

Figure 13(a). Time varying reactance signals at 3.2 mhz for portion of sample with no 

defects. 
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Figure 13(b). Frequency Spectrum 



Figure 14(a). Time varying reactance data for region above cut tows. 



Figure 14 (b) Frequency Spectrum of the cut-tow region 
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Figure 14 (b) Frequency Spectrum of the cut-tow region 



The relationship between time-varying components associated with the probe 
passing over the filaments was anticipated to very straigh-forward and easy to determine 
merely by matching large amplitude frequency components between the various defect 
regions and the region with no known defects. This methodology would have developed 
along the lines used by Doctor, et al;^ 4 except for the real time feature extraction 
required for the experts system approach. Unfortunately that did not work out too well. A 
lot of time was spect in trying to match up frequency components between the various 
spectra, performing an inverse FFT to obtain a noise free time based signal. In this trial 
and error approach only very small repetitive waveforms were removed. The large 
amplitude signals, which make the eddy current signal interpretation difficult, have 
persisted. Therefore it is still difficult to develop a set of rules which will enable an 
expert system to perform eddy current interpretation of graphite epoxy filament wound 
components. 

Another useful concept which the modelling software allows, is to determine 
optimal geometries for the various probes. For example, a series of simulations 
illustrating how geometry affects the conductance measurement of the graphite epoxy 
filament material using the horseshoe probe are given next. 
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Figure 15(a). Finite element model showing longer coil on legs of the horseshoe. 



Figure 15(b). Finite element model showing larger horseshoe(approx. 1.5 times). 



For the cases tried with the same defect geometry, the impedance measurements 
obtained are given in the following table: 
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Case 

Frequency 

Rmax (ohms) 

Rmin(ohms) 

Xmax(o hm s) X mjn (ohms) 

1. Normal 

3.2 Mhz 

4.242 

4.075 

6.991 6.882 

2. longer coil 

3.2 Mhz 

4.068 

3.982 

7.52 7.61 

3. larger probe 

3.2 Mhz 

3.379 

3.562 

7.29 7.50 

4. larger probe 

1.0 Mhz 

1.499 

1.529 

3.15 3.32 

The case identified as Normal is taken from the data shown in Figures 12(a) and (b). 
Note that with the same defect size in layer 2, the following changes in impedance are 

calculated: 





Case 


AR 


AX 

1 


0.167 


0.109 

2 


0.0086 


-0.09 

3 


-0.183 


-0.21 

4 


-0.0030 


-0.17 


From this data, it is obvious that for the same defect modeled here, the largest signals for 
defect identification are obtained from the probe with the larger footprint. Case 4 was 
then run to see if a larger probe might respond better to a lower frequency; however, at 
the lower frequency, the magnetic flux density is too low to provide a good signal for 
the flaw. Hence, for graphite epoxy materials, larger probe sizes should provide larger 
signals for flaw detection. 

CONCLUSIONS 

This research has continued to make in-roads into establishing theoretical 
foundations for the robotic eddy current inspection system in Building 4702 at Marshall 
Space Flight Center. The largest contributor to understanding how eddy current signals 
show respond to known types of defects is still the finite element software. In these 
studies, the finite element software used was "Maxwell", which is a product of Ansoft, 
Inc. 

There is still some difficulty in scanning filament wound materials due to the surface 
roughness of the fabricated parts used in this study. An attempt was made to use FFT 
analysis to pick out signals from these surface effects; however, we were not successful 
in actually performing a good signal clean-up using FFT’s as a post-processing tool. A 
better approach for experiments such as these is to consider using the new Digital Signal 
Processing units being offered in today's marketplace. The goal to develop a stand alone 
expert system for defect analysis will require such technology. Considerations of neural 
nets for cleaning up the noisy signals may also be considered. 
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APPENDIX 


The following program is used to calculate a 3D response from the E field data obtained 
from Ansoft Maxwell models. 

Program RECIPRO.C 

# include <stdio.h> 

# include <math.h> 

# include <string.h> 

# include <dos.h> 

float y [200] ; 
float ansi; 
float ansv; 
int xxll; 
int xxhl; 
int il; 
int numseg; 
float sp; 
float i; 
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void IntegrateVector (float y[], float samper, int xl,int 
xh, float *ans ) ; 

void ClrScr ( ) ; 

void main ( ) 

{ 

int j = 0; 
float ef [200] ; 
float df [200] ; 
float cond; 
float cur; 
float delx; 
float dely; 
float delz; 
char ext [ ] = ".dat"; 
char dfield [200 ] ; 
char efield[200]; 
char dname [20]; 
char ename [20 ] ; 

FILE *dstream; 

FILE *estream; 

ClrScr (); /*clear screen*/ 

printfC Please enter 'E-field without flaw' file name 
: " ) ; 

gets (ename) ; 
strcat (ename, ext ) ; 

printfC Please enter 'E-field with flaw' file name : 

") ; 

gets (dname) ; 

printfC Please enter 'conductivity' of specimen : "); 
scanf("%e", &cond) ; 

printfC Please enter 'source current' : "); 
scanf("%f", &cur) ; 

printfC' Please enter 'defect depth' : "); 
scanf("%f", &dely) ; 

printfC Please enter 'defect width' : "); 
scanf("%f", &delx) ; 

printfC' Please enter 'defect length': "); 
scanf("%f", &delz); 
strcat (dname, ext) ; 

estream = fopen (ename, "r"); 
dstream = fopen (dname, "r"); 


/*read file*/ 

while ( fgets (efield, 64, estream) && fgets (dfield, 64, 
dstream) !=NULL) 

{ 

ef[j] = atof (ef ield+44 ) ; 
df[j] = atof (df ield+4 4 ) ; 
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y[j] = fabs(ef[j] * df [ j ] ) / (cond*cur*cur ) ; 

j = j + 1/ 

} 

/* calculate the impedance change */ 

xxll = 0.0; 

xxhl = 9 9.0; 

sp = dely / 100.0; 

IntegrateVector (y, sp, 0.0, 99.0, &ansl); 
ansv = ansi * delx * delz; 

printf("The Impedance Change = %e\n", ansv); 


/* 1/3 and 3/8 Simpson rule */ 

void IntegrateVector ( float y[], float samper, int xl,int 
xh, float *ans ) 

{ 

int numseg; 
int i; 

int strtpnt; 
int endpnt; 
int segcntr; 
float areal; 
float area2; 
float area; 
char even; 
float segwidth; 

numseg = xh - xl; 
segwidth = samper; 
endpnt = xh; 
areal = 0.0; 
area2 = 0.0; 
area = 0.0; 

if ( fmod (numseg, 2 ) != 0 ) 

{ 

areal = 3.0 / 8.0 * segwidth * (y [endpnt - 3] + 3.0 * 
y [endpnt - 2] + 3.0 * y [endpnt - 1] + 

y [endpnt] ) ; 

endpnt = endpnt - 3; 

} 

else 

{ 

areal = 0.0; 

} 

if ( numseg != 3 ) 

{ 

strtpnt = xl; 
do 
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